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Abstract. Calculations of one-electron spectral functions, optical conductivity 
and spin-wave energy in the Holstein double-exchange model are made using 
the many-body coherent potential approximation. Satisfactory agreement is 
obtained with angle-resolved photoemission results on Lai.2Sri.gMn207 and 
optical measurements on Ndo.7Sro.3Mn03. A pseudogap in the one-electron 
spectrum at the Fermi level plays an important role in both systems, but a small- 
polaron band is only predicted to exist in the La system. A rigorous upper bound 
on spin-wave energies at T = is derived. The spin-wave stiffness constant 
D decreases with increasing electron-phonon coupling g in a similar way to the 
Curie temperature Tq, but D/(k^Tc) increases for large g (low Tq) as observed 
experimentally. 



PACS numbers: 71.27.+a, 71.38.-k, 74.25. Gz, 75.30.DS, 75.30.Vn 
1. Introduction 

Recently there has been much interest in the manganite compounds R 1 _ x A x Mn0 3 , 
where R is a rare-earth element such as La or Nd and A is a divalent metal ion such 
as Ca, Sr or Pb. These compounds are generally ferromagnetic for x ~ 0.2 — 0.4 
and in many of them, near the Curie temperature T c , the electrical resistivity p 
decreases strongly in an applied magnetic field. This effect is known as colossal 
magnetoresistance (CMR) and a recent review is by Ramirez Q. Ferromagnetism 
in these compounds is believed to be due to the "double exchange" (DE) mechanism, 
which operates when local spins are strongly coupled, by Hund's rule, to the spins of 
itinerant electrons occupying a narrow band. The local spins, of magnitude 5 = 3/2, 
correspond to three localized Mn d electrons of t 2g symmetry and the band is derived 
from Mn d states of e g symmetry. The band contains n = 1 — x electrons per atom. 
Millis et al 0] stressed that to describe the physics of the manganites completely 
it is necessary to consider coupling of the electrons to local phonons as well as 
to local spins. Subsequently, Millis et al || made detailed calculations of such a 
model, with local spins and local moments treated classically, using dynamical mean 
field theory (DMFT). Recently Green extending a many-body coherent potential 
approximation (CPA) developed by Edwards et al || ||] for the pure DE model, 
considered a similar model with local spins and local phonons treated quantum- 
mechanically. In the classical spin limit of the DE model the many-body CPA is 
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equivalent to the DMFT of Furukawa |7| || . Green Q discussed a number of physical 
properties of the model with electron-phonon coupling, which he called the Holstein- 
DE model. These include the resistivity, and its dependence on the applied magnetic 
field (the CMR effect) and on pressure, and the Curie temperature T c . In section [2] we 
summarize the many-body CPA treatment of the Holstein-DE model and comment 
on the application of Green's results to the manganites. The aim of this paper is to 
calculate some spectral properties of the model and to compare with experimental 
results. We restrict our attention to the paramagnetic state and the ferromagnetic 
ground state where the many-body CPA simplifies considerably. In section || we 
calculate spectral functions which relate to the results of angle-resolved photoemission 
(ARPES) and in section | we consider the optical conductivity. In section we 
calculate the spin-wave stiffness constant in the ferromagnetic state at T = and 
show how it decreases with increasing electron-phonon coupling in a similar way to 
T c . A brief summary is given in section 



2. Many-body CPA for the Holstein-DE model 

The Hamiltonian of the Holstein-DE model in the absence of an applied magnetic field 
is 

H = E vU* - J E s i ■ °i -»E n i ( b t + 6 +^E fo *V 

ija i i i 

The first two terms constitute the DE model and the first, third and fourth terms 
form the Holstein model ||. Einstein phonons on site i, with energy u> and creation 
operator b\ couple to the electron occupation number n i = ^2 (J n i with coupling 

strength g. Here n ia = c\ a c ia , where c\ a creates an electron of spin a on lattice site i. 
The conduction electron spin a i = (af, of , cr| ) couples to the local spin S i with Hund 
exchange parameter J > 0, and t i - is the band hopping integral. The components of 
er • arc defined as 



H n n- n u)- ( 2 ) 



It is reasonable to assume that the manganites are in the double exchange limit J > W 
Jio| |, where 2W is the width of the itinerant electron band, and in this paper we assume 
J = 00. This ensures that there is no double occupation of a site by electrons, so that 
the system is a Mott insulator for n = 1. This is not so for the two-band model used 
by Millis et al B, and Held and Vollhardt [[ll) have stressed the need to introduce a 
strong on-site Coulomb interaction in this case. Clearly the one-band model neglects 
effects of e g orbital degeneracy such as a proper treatment of the Jahn- Teller effect 
and the orbital ordering which occurs in the undoped system (n = 1). However the 
nature of the local phonon mode in equation (|l|) is not specified and could correspond 
to a tetrahedral distortion of the oxygen octahedron surrounding a Mn site as in the 
dynamical Jahn- Teller effect. 

An important feature of the many-body CPA for the one-particle retarded Green 
function is that it becomes exact in the atomic limit t- = 0. In this limit, and with 
J = 00, the Green function is given by H 



AL _ ~ I r {2A[6H(6H + l)] 1/2 } 
(t " {<i ~ (2S'+l)exp{A[26(w) + l]} 



r — — 00 
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(2S + l)f cxp (r/fa/2) + (5 + 1)(1 - n) exp (-r(3u/2) 

x ; 3 (3) 

e + ru 

where L, is the modihed Bessel function, A = g 2 /uj 2 and b(u>) = (exp(/3a>) — is the 
Bose function with j3 = (fc B T) _1 . In taking the limit J — > oo we have made a shift 
of energy origin e — ► e — JS/2 and the polaron binding energy Aw is also absorbed in 
the chemical potential. The form of the many-body CPA used by Green Q for finite 
band- width becomes particularly simple in the paramagnetic state. For a band density 
of states of elliptic form D (e) — [2/ (nW 2 )}\/W 2 — e 2 , the local Green function G(e) 
satisfies the CPA equation 

G(e) = G Ah (e-W 2 G/4). (4) 

Furthermore the self-energy S(e) is related to the local Green function by the equation 

1 

E(e) = e-G- 1 -W 2 G/A. (5) 

As well as the paramagnetic state, we shall also consider the case of a completely 
saturated ferromagnetic state at zero temperature, with all local and itinerant spins 
aligned. Then the double-exchange term in equation (Q) becomes merely a constant 
shift in energy and the Hamiltonian is equivalent to that of the pure Holstein model. 
Within the many-body CPA the saturated state is actually the self-consistent ground 
state only for S = oo Q and we shall mainly consider this limit. It is found that in the 
DE model JB], |(| neither T c nor the resistivity vary enormously with S so that this is 
a reasonable approximation to the S = 3/2 Mn spin. In the saturated ferromagnetic 
state at T = 0, with all n electrons per atom having f spin, the local Green function 
G| (e) again satisfies equation (^) with G AL replaced by 

Gt L (<) =*- X I 1 +11^7 (— + —)) ■ (6) 
T w 1 e ^-J rl \e + tor e - ujr J J w 

The many body CPA is successful in describing the crossover from weak electron- 
phonon coupling, through intermediate coupling where small-polaron bands begin to 
appear, to strong coupling where some results similar to those of standard small- 
polaron theory are recovered. It therefore extends the work of Millis et al ||, 
where phonons are treated classically, to the quantum small-polaron regime. One 
important quantum effect on a thermodynamic property is the behaviour of T c for 
strong electron-phonon coupling. Millis et al [|| find T c ~ <?~ 4 whereas Green j| finds 
T c is exponentially small for g/W > 0.35. The physics is dominated by a very narrow 
polaron band. Green Q showed how the Holstein-DE model could describe the very 
different behaviour of L ai _ x Sr x Mn0 3 (LSMO) and L ai _ x Ca x Mn0 3 (LCMO), with 
x ~ 0.33 where T c is largest. In LSMO the resistivity p increases monotonically with 
temperature T, with a metal-poor metal transition at T c , whereas LCMO shows a 
metal-insulator transition with p decreasing with T above T c . Also p(T c ) is an order of 
magnitude smaller in LSMO than in LCMO, This type of behaviour can be understood 
in the Holstein-DE model by assuming that g/W is about 50% larger in LCMO than 
in LSMO and, within the model, this is consistent with the considerably lower T c in 
LCMO. (The form of p(T) in figure 7 of reference 4 is indistinguishable for g/W = 0.1 
from that shown for g/W — 0.01 |L2| and contrasts strongly with the curve in figure 6 
for g/W = 0.16). The calculated resistivity is not as sensitive to applied magnetic 
field and pressure as observed but this discrepancy might be removed by introducing 
a decrease in g/W with decreasing resistivity owing to enhanced screening of the ionic 
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charges Q . A defect of the CPA treatment is its failure to describe coherent Bloch-like 
states in the saturated ferromagnetic state at T = 0, leading to a spurious residual 
resistivity Q. 

In this paper we calculate further properties of the Holstein-DE model and 
confront them with experiment. 



3. Angle-resolved photoemission spectroscopy (ARPES) 

The decreasing resistivity of LCMO as the temperature increases above T c is, 
according to Millis et al [|J and Green (J], due to the gradual filling of a pseudogap 
in the density of states. In Green's work the pseudogap contains well-defined polaron 
sub-bands in the (hypothetical) paramagnetic state at T = but above T c , with 
parameters appropriate to LCMO, these are smeared out so as to resemble the classical 
picture of Millis et al . The pseudogap should be observable in ARPES measurements 
and in optical conductivity. In an excellent paper on ARPES for the bilayer manganitc 
La-,^ 2 Sr 1 g Mn 2 7 , nominally with n = 0.6, Dessau et al [ jl3| interpret their results very 
much in the spirit of the Holstein model. In this section we present some calculations 
of the spectral functions in the low temperature ferromagnetic state to compare with 
this experimental data. For convenience we take n = 0.5 and also S = J = oo, as 
discussed above. The charge ordering which might occur for n = 0.5 is suppressed 
since our treatment imposes spatial homogeneity. With this constraint there is no 
qualitative difference between systems with n — 0.5 and with n — 0.3 or 0.6, for 
example. We take the phonon energy to be uj/W = 0.05, the same typical value used 
in previous calculations B. A half- width W = 1 eV is consistent with the e bands 



crossing the Fermi level in the calculations of Dessau et al |13j , shown in figure |3(b) 
The low T c — 126 K in this bilayer manganitc is partly due to quasi-two dimensional 
fluctuations, but the large resistivity p ~ 3 milcm at low temperatures indicates that 
small-polaron bands might exist even in the ferromagnetic state. Consequently the 
clcctron-phonon coupling should be stronger than in cubic manganites like LCMO and 
we choose g/W = 0.2. 

The one-electron spectral function is given by 

A h (e) = -K- 1 lm [e-e^ECe)]- 1 

= - 7r -iS' e 7[( e -e fe -E' e ) 2 + sf] (7) 

where E^, E" are the real and imaginary parts of the self-energy E(e) and e k is the 
band energy for wave- vector k. In the ferromagnetic state at T — the local Green 
function GJe) is calculated from equation (||), with G AL defined by (||), and E(e) 
follows from equation These equations assume an elliptic density of states which 
we here regard as an approximation to the density of states for a band which takes 
the form e k = — W cosiry for k — n(l,y), < y < 1. This band is shown as a full 



line in figure 3(a) and crosses the Ferm i level E F at k = 7r(l, |). It roughly models 



one of the x — y bands in figure 3(b) . The calculated results for A k are shown in 
figure |. Well away from the Fermi level, a well-defined peak exists which broadens as 
k approaches the Fermi momentum at y = 0.5. For larger y the weight below the Fermi 
level is strongly reduced. The peaks never approach the Fermi level closely which is an 
important feature of the observed spectra |l3| reproduced in figure ||. The theoretical 
curves in figure [l] resemble quite closely the data of figure ||(c) . There is a pseudogap 
in the calculated spectra extending about 0.1 eV on each side of the Fermi level. In 
fact this pseudogap contains polaron bands like those shown in figure 4 of Green's 
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paper Q|. However, their amplitude is too small to show up in figure [l] and in the 
experimental data. Nevertheless, it is the central polaron band around the Fermi level 
which is responsible for the low but finite conductivity of the system. The positions 
of the peaks in figure |l| are plotted in figure 3(a) and comparison can be made with 
the right half of figure |3(b)| reproduced from Dessau et al |L3| . Filled and unfilled 
symbols correspond to high and low weights, respectively, obtained by integration 
of the spectral function up to the Fermi energy. This comparison between theory 
and experiment supports the conclusion of Dessau et al Jl3| that, in the manganites 
with a layered structure, strong electron-phonon coupling (with the appearance of a 
pseudogap) is already important below T c . This contrasts with the usual pseudocubic 
manganites where the pseudogap only appears above T c . Previous work related to 
ours is the calculation by Perebeinos and Allen JlJ] of ARPES spectra in a two-band 
model of undoped LaMn0 3 . It should be mentioned that Moreo et al Jll| interpret 
the observed pseudogap not as an intrinsic property but in terms of phase separation. 




Energy relative to E p in eV 



Figure 1. The spectral function A^(t) in the ferromagnetic state at T = for 
J = S = oo, n = 0.5 and strong electron-phonon coupling g/W = 0.20, with 

= f(l:3/)- 



4. Optical conductivity 

In the local approximation of CPA or DMFT there is no vertex correction in the 
current-current response function which may thus be expressed in terms of the one- 
particle spectral function. In the paramagnetic state the optical conductivity is given 

by H | 

k J 

where the electron velocity v k — Ve fc , a 3 is the volume of the unit cell, /(e) = 
{exp[/3(e — fij] + is the Fermi function and N is the number of lattice sites. 

Since A k depends on k only through e k we may define a function 4>{e k ) such that 
4>'{e k ) = A k {t)A k {e + v), where the dependence of on e and v has been suppressed 
in the notation. Hence the sum in equation (|^) may be written as 

E Ve fc- v ^) = -E<K e fc) v2e fc' (9) 

k k 
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Figure 2. ARPES spectra of Lai. 2 Sri. 8 Mn 2 07 (T c = 126 K) in the 
ferromagnetic state at T = 10 K, reproduced from Dessau et al |13|. 
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Figure 3. (a) the band energy ej. ( ) and the position of the peak centres in 

figure |j| Different symbols denote high (a) and low (A) spectral weight, obtained 
by integration over the spectral function up to the Fermi energy Bp. (b) virtual 
crystal LSDA+U jriajority spin hand s fo r Lai.2Sri.sMn207 with experimental 
peaks from figure El reproduced from [n3| . 



the last step following by Gauss's theorem. For a simple cubic tight binding band 
e fc = — 2t J2aS 0S ^a a > w ith a summed over x, y, z, V 2 e^= — a 2 e fe . Then the summand 
in equation (0) is a function of e k only and equation m) becomes 

-(") = ^ / de / 6ED e WEm m -* + V) (10) 



3aS 

where D c (e) is the density of states for the simple cubic band. By using equation (^) 
in the definition of <fi'(e k ), and integrating with respect to e fc , we obtain 

1 1 



£^ 2 + o 2 )((^' + e» + „) 2 + o 2 



K + AK + u-K 2 + 0^)arctanP 
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K 2 - K+v + ° 2 ) arctanQ + E^OtogiZ 



(11) 



with 

= E'-E' +u, Q 



E' 



^' E" 2 +(E' 



Since we calculate the self-energy E using the elliptic density of states D 
discussed in section 0, it is reasonable to approximate D C (E) in equation (J 
D C (E). The integral over J5 can then be carried out by parts so that, using the 
definition of <fi' , we find 

„ M dB^D.C WM.(« + ,) Mz£l±i> (13) 

with A s (e) given by equation (0), e fc being replaced by E. This is of the form given by 
Chung and Freericks jl7[ and Chattopadhyay et al fjsjj . If D C (E) in equation ([[o|) is 
replaced by a Gaussian, corresponding to a hypercubic lattice in infinite dimensions, 
the factor (W 2 — E 2 )/3 in equation ( |l3| ) is replaced by a constant In this 

paper we use equation ([l3]) which is consistent with our previous calculations of 
f7 (0) dj^' This expression satisfies the correct one-band sum rule |l8|, that 
(2/tt) J q cr(^)d^ = — Ke 2 / (3ah), where the "kinetic energy" K is the thermal average 
per lattice site of the first term in the Hamiltonian (|l|) . 

As in section || we take J = S = oo, W = 1 eV, u> — 50 meV and, as previously 
|Q, ||, a = 5A, which is slightly larger than the Mn-Mn distance in perovskite 
manganites. It is again convenient to take n = 0.5 so that for J — S = oo the 
chemical potential fi is fixed at the centre of the occupied band for all temperatures, 
by symmetry. In the ferromagnetic state at T = the spin-degeneracy factor 2 in 
equation ([l^) is omitted. 

There have been numerous optical investigations of the manganites and the 
experimental data and their interpretation are varied. For photon energy v > 3 eV 
various peaks in the optical conductivity have been assigned to interband transitions 
between the Hund's rule split bands and to charge-transfer transitions between the 
O 2p and the Mn e g bands. The latter do not appear in our one-band model and 
our assumption of J = oo eliminates the upper Hund's rule band. Our calculated 
o-(y) is therefore only non-zero in the region v < 2.5 eV and in general exhibits one 
peak. We shall focus the discussion by considering one material, Nd Q 7 Sr Q 3 Mn0 3 
(NSMO), which has been investigated by at least two experimental groups. In the 
paramagnetic state above T r there is one feature common to data on both thin films 
pl[ |20|| and single crystals |pffi| . This is a broad peak at about 1.2 eV with a maximum 
conductivity cr max ~ 0.7 — 0.9 (m£7cm) -1 . However in NSMO films T c , as deduced 
from the maximum in the resistivity p(T), is larger in oxygen-annealed films (T c « 230 
K) than in an unannealed sample (T c w 180 K) |2(], |l|. In a NSMO single crystal, on 
the other hand, p{T) is essentially the same (T c w 200 K) in polished and annealed 
samples but the peak in the polished sample is shifted to v ~ 1 eV with a reduced 
°max ~ 0-3 (mficm) -1 . Since p(T) and a{v) are quantitatively quite similar in NSMO 
and LCMO [^0) we model NSMO with electron-phonon coupling g/W — 0.16, the same 
value as proposed by Green || for LCMO. The dc conductivity cr(0) at T = 10 K is also 
very similar in annealed NSMO and LCMO films, 2.9 and 3.3 (mSlcm) -1 respectively 
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[ pO[ . The value for NSMO agrees well with a measured cr(0) of 3.2 (mficm) 1 at 
T = 15 K in a NSMO single crystal 

Using the parameters discussed above the Curie temperature of the model is 
about 230 K [Q and in figure 4(a) we plot the calculated optical conductivity for the 
ferromagnetic state at T = and the paramagnetic state at T = T c and T = 1.5T C . 
The low dc conductivity of 0.11 (mficm) -1 at T = is due to the inadequate treatment 
of the coherent ground state in the CPA. Fortuitously the incoherent scattering 
introduced by the CPA seems to model quite well the low temperature incoherent 
scattering in the unannealed NSMO sample of Kaplan et al |2l|] which has a dc 
condu ctivit y of about 0.15 (mficm) -1 at T = 15 K. This optical data is reproduced in 
figure \j(b j for comparison with the calculated results of figure 4a. In annealed NSMO 
films |2Cf|, and in single crystals (22), o{y) in the low temperature ferromagnetic state 
continues to rise with decreasing v down to much lower photon energy, and c(0) ~ 3 
(mficm) -1 . The above-mentioned defect of the present CPA treatment is associated 
with the absence of the sharp quasi-particle peak in the spectral function whith should 
exist in the ferromagnetic state at T = 0. Consequently the Drude peak in o{v) at 
low frequency is also absent. 




Photon Energy (eV) 



(a) (b) 

Figure 4. (a) Calculated optical conductivity for strong electron-phonon coupling 

g/W = 0.16 in the ferromagnetic state at T = ( ), the paramagnetic state 

at T = Tq ( ) and the paramagnetic state at T = 1.5Tc ( ). The plot is 

for J = S = oo, n = 0.5 and a = 5A. (b) optical conductivity of Ndo.7Srn.3Mn03 
at different temperatures, reproduced from Kaplan et al fell . 

As pointed out above, cr(v) is much less sample-dependent in the paramagnetic 



state above T c and a quantitative comparison with theory is meaningful. In figure 4(a) 
cr max « 0.15 (mf2cm) _1 and it must be remembered that this curve is the contribution 



of optical transitions within the e g band only. To compare with the data of figure 4(b) 
one should subtract the background due to other transitions which might leave an 
effective cr max w 0.3 (mf2cm) _1 . Bearing in mind the simplicity of the one-band model, 
this order-of-magnitude agreement between theory and experiment is satisfactory. Our 
calculated results are quite similar to those of figure 7d in Millis et al (^] . Clearly, for 
the present intermediate electron-phonon coupling strength, their classical treatment 
of phonons is sufficient to obtain the essential features of the optical conductivity. 
The calculated peak in a{v) with v k, 1 eV arises from fc-conserving transitions across 
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the type of pseudogap discussed in section During the process an electron moves 
from one site to a neighbouring one which was previously unoccupied. The electron 
motion is accompanied by a lattice distortion, of Jahn- Teller type, which corresponds 
to a displacement of the local phonon oscillator coordinate in the Holstein-DE model. 
When an electron enters (leaves) a site the final displaced (undisplaced) oscillator is 
generally in an excited state with typical excitation energy g 2 /ui. This is the atomic- 
limit polaron binding energy and for the parameters assumed here is about 0.5 eV. 
Thus the peak in cx(z/) occurs at about twice the polaron binding energy just as in the 
standard small-polaron theory (^3). However for the present intermediate electron- 
phonon coupling g/W = 0.16 polaron bands near the Fermi level are largely washed 
out above T c j|, so standard small-polaron theory is not expected to apply to a{v) 
for low photon energies. In fact an activation energy in the dc conductivity of 0.25 eV, 
half the polaron binding energy as predicted by small-polaron theory []23[ , is about 
a factor 4 larger than one deduced from Green's numerical calculations. Green's 
calculation of the dc resistivity above T c is in reasonable agreement with experiments 
on NSMO where an activation energy of about 0.08 eV is found Lee et al Q 
quote an activation energy of about 0.15 eV and, with undue reliance on small-polaron 
theory, expect to have a peak in a{v) at about 0.6 eV. Although there is no sign of 
such a peak in their data they claim that their one broad peak near 1.2 eV should be 
interpreted in terms of a two-peak structure, one near 1.5 eV and the other below 1 
eV. Our interpretation of the 1.2 eV peak in a(y) for T > T c is broadly in line with 
that of several other previous authors ||, ^(| |21[] . However we stress again that, for 
moderate clcctron-phonon coupling standard small-polaron theory does not hold at 
low photon energy, where only states around the Fermi level are involved, so no simple 
link between peak position and the activation energy of the dc conductivity can be 
made. Below T c the peak in aiy) shifts to lower frequency as the pseudogap rapidly 
fills in. However in our calculations this shift is held up due to spurious incoherent 
scattering in the ground state, which limits the low-temperature dc conductivity. The 
same effect actually occurs in unannealed NSMO films (figure 0) but in annealed films 
and single crystals the peak shifts almost to zero frequency p2[ |. The theoretical 
situation could be improved by introducing screening effects in the electron-phonon 
interaction so that g/W decreases with decreasing resistivity, as mentioned at the end 
of section ||. The spurious residual resistivity at T — drops from 9 mficm to 0.5 
mficm as g/W decreases from 0.16 to 0.10 J|, [l^]. Ishihara et al and Mack and 
Horsch |26|] have given an alternative interpretation of the broad low energy peak at 
T = in terms of orbital degrees of freedom in the doubly-degenerate e g band. They 
propose that strong correlation, with the constraint of no doubly occupied sites, leads 
to incoherent motion of the carriers. 



5. The spin- wave spectrum 

Quijada et al |Q considered the spin- wave stiffness constant D in the saturated 
ferromagnetic state at T = and its relation to optical conductivity. However their 
derivation of an expression for D was restricted to the DE model. Here we derive 
a simple approximate formula for the spin-wave spectrum of the Holstcin-DE model. 
Our main aim is to calculate D as a function of clcctron-phonon coupling strength 
and to compare with the corresponding behaviour of T c . In the limit of J — > oo the 
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local and itinerant spins are locked together as we make the ansatz 

\Q)={s-+a-)\F) (14) 

for the state (unnormalized as yet) with a magnon of wave-vector q excited. Here \F) 
is the exact ferromagnetic ground state, which we assume to be one of complete spin 
alignment, and the spin lowering operators are defined by 

s , = E eiqRiS i > % = E ei9i V > ( 15 ) 

i i 

where i?^ is the position of lattice site i. 

Actually we do not know the state \F) exactly for the Holstein-DE model, unlike 
the DE model, but we shall only approximate it at a later stage. The spin- wave energy 
is given by 

(gjgjg) 
(q\q) 

where E Q is the exact ground state energy, so that H \F) = EAF) . By the variational 
principle, this is an upper bound on uj q . Using equation ( |l4|) it is easy to show that 

(q\[H,S-+a-]\F) 

w « = — re — (17) 

and (q \ q) — N(2S + n), where 5 is the magnitude of the localized spin S, r Only the 
first term in the Hamiltonian (Q) makes a non-zero contribution to the commutator. 
This one-electron term may be written as ^2 ka £ k n kcJ where n ka is the occupation 
number for the Bloch state ka. A straight-forward calculation yields 

= N(n + 25) £ ( £fc+ " " £fe ) (18) 

where (n k ^) = (F\n k ^ \F). Assuming a simple cubic tight-binding band e k = 
—t^2 R e lk R , the sum being over 6 nearest neighbour sites, we find 

Since (n k ^) has cubic symmetry in fc-space, the k sum in equation Jl9| ) is independent 
of the particular neighbour R. By including the factor —t it may therefore be written 
as (l/6)£ fe e fe (n feT ). Hence 

a R 

— (2D /a 2 ) ^3 — cosg^a — cos q y a — cos q z aj , (20) 
where q= (q x ,q y ,q z ) and 

a 2 v ^ Ka 2 

U = ~ 6N(n + 25) ^ €k{n ^' = ~ 6(n + 25) ' ( ] 

Here K is the expectation value of the kinetic energy which appears in the optical 
sum rule mentioned in section |]. By expanding equation (|2C| ) in powers of q to second 
order, we find ui — Dq 2 as that D defined by equation ( |2lD is the spin-wave stiffness 
constant. 
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At this stage, with the expectation value (n k «) calculated in the exact ground 
state, equation ( pp[ ) is a rigorous upper bound on the magnon energy for arbitrary J 
and S. This is no longer the case when we proceed to evaluate it within the many-body 
CPA. Clearly, equations ( |2C| ) and ( pl| ) apply equally to the Holstein-DE and the DE 
model. However, we shall show that electron-phonon coupling in the former model has 
a strong influence via (n^). In the limit J — ► oo, equation ( pl| ) is equivalent to results 
derived by Nagaev [§7[ §§], Kubo and Ohata |§, Furukawa §, |o|, Wang §l| and 
Quijada ei aZ [^(j for the DE model. The K defined in the latter paper is one-third 
of our kinetic energy. Furukawa also derived the dispersion relation of equation ( pp|) 
for the DE model. Now 

k t > = r deA ^) (22) 

J —oo 

where A k (e), given by equation (]?]), is calculated for the saturated ferromagnetic state 
at T = 0. As in deriving o~(v) we can replace the k sum in equatio n (|2l| ) by an energy 
integral. Hence, using the notation A E (e) introduced in equation (|13|), 

D = L AE L •" ED -^ A ^- < 23 » 

As in the previous section we have approximated the simple cubic density of states by 
the elliptic density of states D e (E). To match the other calculations in this paper, in 
particular that of the optical conductivity, we calculate A E (e), and hence the double 
integral in equation (|2^) which represents the average kinetic energy K, in the limit 
S = oo. However, in the prefactor in equation (^3|) we put S = 3/2 as is appropriate 
for the localized Mn spins. This corresponds to calculating the Bloch wall stiffness 
constant (oc D(n + 25)), which is a static quantity, in the limit S — > oo but retaining 
the essential finiteness of the spin in the dynamical quantity D. Otherwise we use 
the same parameters as in previous sections, except that we now take a — 4 A. In 
section |] we used o = 5A, for consistency with our earlier work on conductivity, but 
4 A is closer to the Mn-Mn distance in the pseudocubic manganites. 

In figure || we plot the spin-wave stiffness D at T = as a function of 
electron-phonon coupling g/W. The reason for the striking decrease of D with g/W, 
particularly in the range 0.1 < g/W < 0.2 applicable to the manganites, is clear 
from equation j2l|). For g/W = 0, the pure DE model, (n k *) = 1 for k within the 
Fermi surface and (n k ^) = otherwise. The negative quantity K is the full non- 
interacting one-electron energy of the ferromagnetic state which drives the double- 
exchange mechanism. For larger g/W, (n k ^) is more spread out over the whole zone 
and \K\ decreases. In an extreme limit where electrons are localized at sites, (n k ^) 
is constant throughout the zone and hence D = 0. This behaviour of D in the 
Holstein-DE model is very similar to that of T c , as calculated by Green pj]. The main 
difference is in the extreme strong-coupling limit where T c becomes very small at 
g/W w 0.35 whereas D is decreasing quite slowly. The slow decrease of D is exactly 
what one expects from equation ( |2l| ) and small-polaron theory, where the kinetic 
energy K ~ g~~ 2 J32|. T c seems to be determined more by the width of the narrow 
polaron band around the Fermi level, which decreases exponentially with g. Thus one 
may expect that D/(k B T c ) increases with increasing g, and thus with decreasing T c . 
This is found experimentally, as discussed later. 

Green and Edwards |(| find that, in the DE model with n = 0.5, T c only increases 
by 5% when S increases from 3/2 to oo. A similar insensitivity to S for 5 > 3/2 is 
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Figure 5. The spin-wave stiffness D versus electron-phonon coupling g in the 
saturated ferromagnetic state at T = 0. The plot is for S = J = oo, W = 1 eV, 
n = 0.5 and a = 4 A. 



expected in the Holstein-DE model. Hence, from equation (23), with A E (e) taken 
in the 5* — » oo limit as discussed, we have 6 = D / (k B T c a 2 ) oc (S + in) -1 . This 
is similar to the result 6 = ^(s + 1) _1 for the spin s nearest neighbour simple cubic 
Heisenberg model, with T c calculated in mean field theory. Using values of T c accurate 
to within about 1% we find an improved Heisenberg value of 5 = 0.286 for s = 3/2, 
and 8 — 0.258 for an interpolated s — 1.75 which models the spin per site S + \n 
in the present model. In the present calculations for D/a 2 , together with those of 
Green (§ for T c , we find, for S = 3/2, S « 0.24 for g/W = 0.1 and 6 w 0.29 for 
5 /W = 0.16. With a = 4 A these values correspond to D/(k B T c ) m 3.9 and 4.6 A 2 , 
respectively. The agreement of these values of 5 with the Heisenberg model suggests 
that for moderate g/W the ferromagnetic transition in the many-body CPA treatment 
of the Holstein-DE model is quite Heisenberg-like. However, this may not be the case 
for larger g/W where, as discussed above, 5 increases rapidly with g/W. 

To compare our theoretical results with experiment we first note that the simple 
spin- wave dispersion in equation (pp|), which is of the Heisenberg form, has been 
found to fit data on La 7 Pb 3 Mn0 3 (T c = 355 K) throughout the Brillouin zone 
[ ^5| . In this work the low temperature stiffness constant D = 134 meVA 2 so that 
D/(k B T c ) = 4.4 A 2 . In LSMO, also with x = 0.3, Martin et al (33) find T c = 378 K 
and D » 188 meVA 2 (at 27 K) so that D/(k B T c ) = 5.8 A 2 . Pr Q 63 Sr 37 Mn0 3 seems 
to behave similarly with T Q = 301 K, D = 165 meVA 2 , D/{k B T c )'= 6.4 A 2 @. 
Although D/(k B T c ) in the last two materials is larger than in a Heisenberg model, 
their spin dynamics near T c is quite conventional |33|, |3(| . However in some systems 
with lower T c this is not the case, and D/(k B T c ) is considerably larger. Thus in 
LCMO, with T c = 250 K, D = 170 meVA 2 , D/(k B T c ) = 7.9 A 2 , the ferromagnetic 
transition seems not to be a standard second order one [[37J. Also in NSMO, with 
T c = 198 K, D = 165 meVA 2 , D/(k^T G ) = 9.7 A 2 , the spin-wave stiffness constant 
does not collapse to zero at T — T c |36|| , just as in LCMO |37|]. It is not clear whether 
such behaviour could occur within the Holstein-DE model or whether inhomogeneity 
due to disorder is important. However the larger values of D/(k B T c ) predicted by 
the model for strong electron-phonon coupling suggests that something unusual is 
going on. It is worth mentioning that D may be underestimated by the many-body 
CPA for an intermediate electron-phonon coupling such as g/W — 0.16 appropriate 
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to LCMO. This is because spurious incoherent scattering near the Fermi level reduces 
the "kinetic energy" K, and hence D. On the other hand equation ( |2l| ) itself gives an 
upper bound to D which, in a better approximation will certainly be reduced. For the 
DE model, in approximations equivalent to the random phase approximation (RPA) 
[ [ju| , |3l| , p0[ , this is achieved by an additional negative term proportional to A -1 where 
A = JS is the Hartree-Fock exchange splitting between up and down spin bands. In 
the presence of an onsite Coulomb interaction U, A — JS + Un so that the negative 
term is considerably reduced. RPA estimates of D in a two-band model differ widely 
po| . For the one-band model with J = oo and S = 3/2 Golosov |59| has shown 
that D is reduced to about half its 5 = oo value, over a wide range of band-filling, 
when magnon-electron scattering processes are considered. Further work on D in the 
Holstein-DE model is highly desirable. 

6. Conclusion 

The many-body CPA treatment of the Holstein-DE model has been used to investigate 
several spectral properties which may be compared with experimental data on the 
manganites. We have been able to supply the theory which was hinted at by Dessau 
et al |fl3f in the discussion of their ARPES measurements on the layered manganite 
La-,^ 2 Sr 1 g Mn 2 7 in the low temperature ferromagnetic state. Broad spectral peaks 
lie either side of a pseudogap at the Fermi level and the pseudogap contains polaron 
subbands with exponentially small weight. One of these, at the Fermi level, is 
responsible for the poor metallic behaviour. We therefore agree with Alexandrov 
and Bratkovsky ]40| that in this system, with unusually strong electron-phonon 
coupling, small polarons exist in the ferromagnetic state. However we find that 
small-polaron theory does not apply above or below T c in a pseudocubic manganite 
like Nd Q 7 Sr Q 3 Mn0 3 with intermediate coupling strength. In particular the small- 
polaron result that the activation energy of the high-temperature dc conductivity is 
one-quarter of the peak photon energy in optical conductivity a(y) is found not to hold, 
in agreement with experiments on NSMO. The observed shift in spectral weight of a(y) 
to lower energy on going into the ferromagnetic state is found to occur, although it 
is somewhat suppressed by spurious incoherent scattering at T — which is a defect 
of the theory. A rigorous upper bound is derived for spin-wave energies at T — in 
the Holstein-DE model. It is shown that the spin-wave stiffness constant D decreases 
with increasing electron-phonon coupling strength in a similar way to T c . However, for 
strong coupling the ratio D/(k B T c ) increases quite rapidly with increasing coupling 
strength, i.e. with decreasing T c . This trend is found experimentally. 
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